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SUMMARY 

The coupled bendlng-bendlng-torslonal equations of dynamic motion of 
rotating, linearly pretwisted blades are derived Including large precone, sec- 
ond degree geometric nonllnearltles and Coriolis effects. The equations are 
solved by the Galerkln method and a linear perturbation technique. Accuracy 
of the present method Is verified by comparisons of predicted frequencies and 
steady state deflections with those from MSC/NASTRAN and from experiments. 
Parametric results are generated to establish where Inclusion of only the sec- 
ond degree geometric nonllnearltles Is adequate. The nonlinear terms causing 
torsional divergence In thin blades are Identified. The effects of Coriolis 
terms and several other structurally nonlinear terms are studied, and their 
relative Importance Is examined. 


INTRODUCTION 

It Is now widely recognized that the Inclusion of geometric nonllnearltles 
In the equations of motion of rotating elastic blades Is necessary for a fair 
prediction of their dynamic characteristics. For the analysis of helicopter 
rotor blades, an appropriate set of nonlinear equations based upon Euler- 
Bernoulll theory was found to yield satisfactory results since the blades were 
essentially slender (refs. 1 to 5). However, there remain certain questions 
concerning the degree to which the geometric nonllnearltles should be retained, 
and concerning the Initial assumptions In prescribing an ordering scheme 
(refs. 6 to 7). Another area, somewhat similar dynamically, but further com- 
plicated due to the geometry, Is the advanced turboprop blade dynamics 
(refs. 8 to 9). In contrast to helicopter rotor blades, the turboprop blades 
are more plate or shell-like, possess variable sweep along the span, and have 
smaller thickness ratios. Furthermore, the ratio of the rotational speed to 
the first nonrotating normal mode frequency of the turboprop blade Is usually 
less than one, while the corresponding ratio for the helicopter blade can be 
somewhat more than one. Because of the unusual geometric features of the 
advanced turboprop blades, finite element methods are normally used. However, 
when finite element modeling Is used. It Is difficult to obtain a physical 
understanding of the complicating effects. 

In order to determine the Individual and combined effects of pretwist, 
sweep and rotation on the blades, parameterlc studies are conducted by using a 
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simpler beam model In which the effects of sweep are Incorporated by preconing 
the blade with respect to the plane of rotation. Precone, a component of 
sweep. Is of the order of 50° for an advanced turboprop blade at the tip. 
Assuming constant precone, linear pretwist and steady rotational speed, the 
coupled nonlinear equations for bendlng-bendlng-torslonal motion are derived 
by using the theory presented In references 3 and 4. In deriving the equa- 
tions, It Is assumed that the elongations and shears are negligible compared 
to unity and the squares of the derivative of the extenslonal deformation of 
the elastic axis Is negligible compared to the square of the bending slopes. 
Shear deflection and rotary Inertia effects are not considered. Geometric non- 
linearities are retained up to second degree. Parametric results from the 
present beam theory for the special case of torslonally rigid blades (refs. 10 
and 11) Indicated excellent agreement with finite element model results 
generated by MSC/NASTRAN. 

In the present paper, the Important effect of torsional flexibility Is 
addressed. A brief derivation of the coupled nonlinear equations of bendlng- 
bendlng-torslonal motion are first presented. Next, Implementation of the 
Galerkln method with nonrotating normal modes for the solution of the nonlinear 
steady state equations and the linearized perturbation equations Is presented. 
Parametric results generated from the present beam theory are compared with 
similar results produced by MSC/NASTRAN. The limitations of restricting the 
geometric nonl Inearltles to second degree, and the effects of Coriolis forces 
for blades with various thickness ratios are shown. Results obtained from 
experimental tests are also presented for typical precone angles, rotational 
speed and setting angles. Finally, the accuracy of the present equations, the 
limits where the second degree geometric nonllnearltles are adequate, and the 
parameters affecting the onset of Instability are discussed. 


EQUATIONS OF MOTION AND METHOD OF SOLUTION 

Figure 1 shows a linearly pretwisted, preconed, and rotating blade of 
uniform rectangular cross section. The coupled bendlng-bendlng-torslonal 
equations of motion for such a blade are derived by using the theory presented 
In reference 3, and by Including large precone and linear pretwist over the 
blade length. When shear and rotary Inertia effects are Ignored, and second 
degree geometric nonllnearltles and Coriolis effects are retained, such equa- 
tions reduce to the following form (a list of notation Is given In appendix B) : 
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Torsion: 
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In deriving the equations given above, an ordering scheme given In ref- 
erence 3 Is followed. Geometric nonlinear terms up through second degree only 
are retained. In general, terms up to 0(e 4 ) In the elastic forces and 0(e 2 ) 

In the Inertial forces are needed In the bending equations. For the sake of 
verification, certain linear terms of the order 0(e 2 ) in the Inertial forces 
are Included In equations (1) and (2). However, these additional terms have 
been found to have an Insignificant affect on the solutions and can be dis- 
carded safely. Considering the torsion equation, terms of the order 0(«5) 

In the elastic forces and terms of the order O(c^) In the Inertial forces 
forces together with the structural and Inertial warping terms, [viz. ( EC-j <t>" ) " 
and (mkj[Q 2 cos 2 0p C <♦>')'], are retained. A large number of higher order terms 
other than those shown In the second degree torsion equation are believed to 
be unimportant (see ref. 3), and are thus discarded. The well known tennis 
racquet effect term appears In the torsion equation as mQ 2 <$> cos 2 Bpcfkjfe - 
kfo) cos 2e, while the tension-torsion coupling term appears In 
[ EAk^u ' ( ©pt + ♦')]'• The extenslonal degree of freedom In the present equa- 
tions has been discarded, since It has been established In reference 11 that 
the effect of the extenslonal coupling on vibration and stability character- 
istics for practical blade configurations Is Insignificant. One can eliminate 
the extenslonal slope, u', by using equations (4) to (6). The effects of the 
tension-torsion and tennis racquet terms and the structural and Inertial 
warping terms were discussed In detail In references 12 and 13 for uncoupled 
torsional vibration of pretwisted, rotating blades. In view of the discussion 
presented In these references, the Inertial warping term Is discarded and the 
structural warping effect terms Is retained. 

Defining the following parameters, 


w = w/L, v = v/L, n = x/L, t = Qt, R = R/L, etc., (8) 

assuming solutions are separable In time and space, and making note of the 
following relations 


d _ dn d _ld d d 

dx " dx dn “ L dn’ dt " Q dx etc -’ 

one can write equations (1) to (3) In the following nondlmenslonal forms: 


(9) 


5 



w + 2 sin 8 v - w sin 0 „ + 
pc £C 


2 cos 6 4 » (mk 2 _ sin 2 e + mk 2 cos 2 e' 
pc \ m2 ml 


f n _ 2 2 

- j sin 0 pc cos 0 pc J ( v 1 + w') dn 


1 



mL 


pc 


- 2 cos B n „ 

r 1 : _ 

w" 1 v dn - w'v 

pc 

J 

L n 


+ sin cos 0_„ 
pc pc 


Jw" j w d n - w'wj 


♦ tw 1v (cos 2 e ♦ sin 2 e 
d 2 


+ w'" (2 y 5 sin 20)^ " 1 j + 2y? * C0S 26 (^2 " 1 

+ ijv - 1 ^ sin 0 cos 0 + v 11 1 (2y5 cos 20) ( - 1 

\d 2 / \d 2 ) 

- v" (2y 2 Sin 20) " 1 ) + 5 (^2 " 1 ) ^ Q % s1n 20 


+ 2w lll 4> 1 sin 20 + w'V sin 20 + 4yw‘ 1 '<)> cos 20 + 4YW ll <t> 1 cos 20 

- 4y 2 w M <t> sin 20 + v cos 20 ± 2v lll (t> 1 cos 20 + v'V 1 cos 20 

~ : r i 

- 4yv 1 1 '<!> sin 20 - 4yv ll <t> 1 sin 20 - 4y v"<t> cos 20 ( 


GJ 


n*A 4 


(4> ' v ' ‘ ' + 2 <t>"v" + 4,'" v ' } 


Y 1 / h 2 \ 

- n sin 0 cos 0 - — SM-r- 1 sin 0 cos 0 sin 20 

pc pc al 2 l rf 2 I pc pc 


( 10 ) 


6 




7 




? : - / 2 

<*> + f, cos 0 (<(> cos 20) + 2f n cos 0 v sin 6 cos 0 + 2 cos 0 w f_ sin 0 

1 pc 1 pc pc \ 2 


+ f 3 cos 2 ©) + f 4 4> 1v ♦ sin 0 pc cos 0 pc 




dn - 4>'w - yvi 


• • 


■ cos V (4,,,Q ” <t> S ^ " 2 C0S ■ V dn - 4>‘V - YV 


+ f ( .(2Y<J> l 4> n + y V) - - f 7 <t> M + f 7 (v‘w ,M + v-w«) 


f f Q cos B„<J>" - 2yf Q sin B nn cos B w" - f, sin B„ cos B v* cos 
8 pc^ 8 pc pc l pc pc 


20 


- - - 2 - 2 

f. sin B cos B„ w' sin 20 + f 0 v"w M cos 20 + f Q (w" - v" ) sin 0 cos 0 
1 pc pc _y yj 


2 2 

f, cos 0 sin 0 cos 0 - yS cos 0„„ 

1 pc 'pc 


( 12 ) 


where 


Q - R(1 - n) + 0.5(1 - n ), S = (R + n), 

, k m2 - k m1 b 2 - d 2 , C b 2 

1 = ,2 ~ i_2 . j2’ r 2 ' , 2 " u 2 . j2* 


m 


b + d* 


m 


b + d 


H 2 EC, 

d - f = 1 

l » 1 A 


2 2 
EbV 


ml 

3 = k 2 - b 2 ♦ d 2 ’ ‘ 4 " m* 2 k 2 L 4 ' 12^ 2 L 4 (b 2 ♦ d 2 ) 
m m ' ' 


. EAk A E(b 2 ± d 2 ) . EB 1 9 it 

f 5 = _2,2,4 " ,4^2 * r 6 " 2. 2. 2 = 


Ey 2 b 4 


mfi k„L 12pL ft 
m 


l | l , l n c 2 . 4 / 1 2 .2! 

k m L 1 5pQ L ( b + d ) 


f GJ _ 4Gt 2 , mk X 

r 7 = . 2 . 2 , 2 “ „ 2 . 2, .2 . J 2 ’ 8 = 2 . 2 


b 2 t 2 


m£2 k^L pft L (b + d") 


mk"L" 1 2L 2 (b 2 + d 2 )’ 
m ' ' 


8 



f„ = 


E(I ^- V 


E(b 2 - d 2 ) 


P Q 2 L 2 d u + I nn ) p« 2 L 2 (b 2 + d fc > 


2 ’ f 10 


GJ 


2 4’ 
msri/ 


El 


11 


5 o-1.5« 


JUL 


• * 4_2 
pAL Q 


, 0 = a + yn, 0_ + = r~ = Yn, 


_ XX _ 
'pt " L “ 


w' = ^ (w) and w = ^ (w) . . 


( 13 ) 


Before discussing the method of solution, It Is worthwhile to point out 
the various Important linear and nonlinear terms In the present equations. In 
equations (10) and (11), certain linear and nonlinear terms are underlined 
once. These terms constitute the linear and nonlinear Coriolis force terms and 
nonlinear terms arising from the foreshortening, (up), and tension terms, 

(Tw')‘ and (Tv')'. Individual and collective effects of these Important terms 
were assessed In references 10 and 11. It was shown In these references that 
the linear and nonlinear Coriolis effect terms significantly affect the onset 
of Instability for large preconed thick blades and thus they must be retained 
In analyzing blades having moderate to large thickness ratios. It was also 
shown that the nonlinear terms In the coupled bending-bending equations which 
vanish for zero precone are extremely Important. These terms can produce sig- 
nificant frequency changes (of the order of 20 percent Increase In the funda- 
mental mode frequency for a 45° preconed thin or thick blade rotating at a 
speed equal to the fundamental frequency of the same blade). Next, the 
Coriolis effect terms due to torsional coupling In the bending equations and 
the flexural coupling In the torsion equation are addressed. Referring to 
equations (10) and (11), one can find the Coriolis effect terms 
[2 cos B pc 4 >(mkjn 2 sin 2 0 + mk^ cos 2 e)/mL 2 ]' and [2 cos & pc <i>(mk(n2 - mk Sl ) 
sin 0 cos 6/mL 2 ]'. The corresponding terms In the torsion equation are 2f q V 
sin 0 cos e cos B pc + 2 cos Bp C w'(f 2 sin 2 0 + f 3 cos 2 0). It may be noted that 
the linear Coriolis forces, 2 sin B pc v and -2 sin Bpc^. the two bend- 
ing equations Indicate skew symmetry of the gyroscopic matrix. One may thus 
consider that Inclusion of torsional degree of freedom does not alter this 
nature of the gyroscopic matrix provided that an appropriate sign change Is 
made In the torsion equation throughout. Keeping In view the ordering scheme 
followed, the bending equations should contain terms of order 0(e 2 ) In the 
Inertial terms. Consequently, the Coriolis effect terms In the bending equa- 
tions which are associated with the torsion variable 4> should be discarded. 

To preserve the skew symmetric nature of the gryoscoplc matrix In linear 
Coriolis force terms, one must discard the corresponding Coriolis force terms 
In the torsion equation also. In the following, we maintain this consistency 
and disregard the aforementioned Coriolis effect terms associated with <j> 

In the bending equations (10) and (11) and all the Coriolis force terms In the 
torsion equation. 


Finally, the twice underlined terms In the bending equations are the non- 
linear terms arising from torsional coupling with the bending motions (pitch- 
flap and pitch-lag couplings). The corresponding terms In the torsion equation 
are shown by underscoring them thrice. These terms were first discussed by Mil 
et al. (ref. 14) and subsequently by others. From the present equations, It 
can be seen that these terms are premultlpl led by the rotational parameter 
E and square of thickness ratio factor, (b 2 /d 2 - 1). For thin blades, the 
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flatwise steady state deflection, w, will become significant with Increasing 
speeds and makes these terms quite significant In the perturbation equations. 
Furthermore, these effects are aggravated due to the coefficient (b 2 /d 2 - 1) 
as the thickness ratio (d/b) becomes smaller. Thus, one may anticipate a 
significant Influence of these terms on the vibration and stability behavior 
of thin blades subjected to large rotational speeds. The terms shown by the 
dashed underlining are associated with the torsional rigidity, GJ, and the 
terms shown by the double dashed underlining arise due to tension effects. It 
will be shown that these terms are not as significant as the aforementioned 
nonlinear terms. 

The coupled bendlng-bendlng-torslon equations are solved by the Galerkln 
method by expanding the dimensionless deflections In terms of a series of gen- 
eralized coordinates and mode shape functions as follows: 


- j < W Oj . *W j <T))„j<l,) 

(14) 

-? (v oJ * 

(15) 

=^(<l> 0 j * A<t>j('0)ej(n) 
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where 

<(/j(n) = cosh (Bjii) - cos ( B^n) - <*j[s1n O^n) - sin (B^n) ] (17) 

6j(n) = 2 sin (Y^n) (18) 

Yj - v(j - \) (19) 

Equations (17) and (18) represent the nonrotating normal modes for a 
cantilevered beam fixed at n = 0 and free at n = 1. The values of Bj 
are taken from reference 15. It may be noted here that sinusoidal mode J 
shape assumed for the torsional degree of freedom Is not compatible with the 
boundary conditions when warping Is Included. However, the effect of warping 
Is not significant for large aspect ratio blades (L/b > 6 (ref. 13)). Thus, 
the mode shape assumed here should produce satisfactory results for blades with 
moderate aspect ratios. The quantities w 0 j , v 0 j and 4 > 0 j are the equilib- 
rium quantities while Awj, Avj and A<t>j are the perturbation quantltltes 
In the generalized coordinates. 

Proceeding as In reference 11, one can apply the Galerkln process for the 
solution of the nonlinear steady state equations and the linearized perturba- 
tion equations (expressed In terms of the equilibrium generalized coordinates) 
which define the blade motion about the equilibrium operating condition. 

The steady state equilibrium equations, and the linear perturbation 
equations are written In the following forms for a solution with n normal 
modes assumed for each Independent variable: 
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with J., LN being the linear and nonlinear parts of the equilibrium equations, 
and JHT^and JK being the mass, Coriolis and stiffness matrices respec- 
tively. ~Element7 of those matrices are presented In appendix A. 


RESULTS AND DISCUSSION 

The nonlinear, steady state, equilibrium equations (20), and the eigen- 
value problem that results from the transformation of equation (21) were 
solved by using computer programs developed In FORTRAN language. Integrations 
are performed on the computer using a 15 point Gaussian quadrature formula. 

The general computer program developed for the solution of equation (21) gives 
the natural frequencies per unit rotational speed, (p/Q). In the presence 
of Coriolis forces, the frequencies will occur In pairs of purely Imaginary 
quantities for a conservative system. In the absence of Coriolis effects, the 
frequency equation (21) reduces to a standard eigenvalue problem the eigen- 
values of which are real quantities, (p^/n 2 ), Specialized cases were solved 
by modifying the general computer program. Parametric studies were conducted 
for thickness ratios (d/b) ranging from 0.05 (representing approximately an 
advanced turboprop blade) to 0.25 (representing approximately a conventional 
propeller blade) for various precone angles, pretwists, setting angles and 
rotational speeds. The disc radius Is assumed to be zero In most of the 
calculations except for cases corresponding to a test configuration developed 
at the NASA Lewis (ref. 16). Results were also generated from the finite 
element code, MSC/NASTRAN, using 250 or 500 CQUAD4 elements for the purpose of 
comparisons with present theoretical results. It should be noted here that 
the MSC/NASTRAN calculations Ignore the Coriolis effects although there Is no 
restriction on the degree of geometric nonlinear terms. All these results are 
presented and discussed In what follows. 


Convergence 

The convergence of the solutions produced by the Galerkln method with 
various numbers of nonrotating normal modes In the coordinate functions Is 
Illustrated In table I. The blade considered for this convergence study has a 
precone of 15°, pretwist of 30° at Its tip, and a thickness ratio of 0.05. 

The blade chord at the root Is set perpendicular to the axis of rotation 
(a = 0°) and the blade rotational speed Is one-half of the fundamental mode 
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frequency of the same nonrotating blade with zero pretwist and zero precone 
(Q/gv| = 0.50). The frequency ratios (p/\) shown In table I are representative 
of the convergence trend that one might expect for a general case of pretwisted 
rotating blade when geometric nonl Inearltles up to second degree and Coriolis 
effects are Included In the analysis. For the purpose of comparison of the 
present beam theory results, and to provide a measure of accuracy of the pre- 
sent beam theory formulation, the results produced by MSC/NASTRAN by using 500 
CQUAD4 elements are also Included In this table. The convergence pattern of 
the components of steady state, dimensionless, tip deflections are shown In 
table II for this blade as obtained from the Galerkln method calculations 
together with those produced by MSC/NASTRAN. 

From the convergence pattern of the frequency ratios presented In table I, 
It can be seen that a five mode solution produces the lowest six coupled 
bendlng-bendlng-torslon frequencies that are In reasonable agreement with the 
corresponding values produced by MSC/NASTRAN. Furthermore, the steady state 
deflections produced by the present beam theory calculations with a five mode 
solution agree quite closely with those from MSC/NASTRAN. Among the several 
factors unknown so far, the sinusoidal mode functions used for torsional 
deformation Instead of hyperbolic functions may perhaps be responsible for the 
oscillatory convergence trend observed for the frequencies and steady state 
deflections shown In tables I and II as the number of nonrotating normal modes 
In the assumed solutions are Increased. It emerges clearly, however, that the 
convergence of higher mode frequencies to accurate values can be accomplished 
by Increasing the number of modes In the solution, and that accuracies of 
practical Interest can be achieved from the present theoretical formulation. 
Finally, the close agreement of the present theoretical results. Including 
Coriolis effects, with the MSC/NASTRAN results, that Ignore the Coriolis 
effects, substantiate the conclusion reached by the present authors In 
reference 11 that Coriolis effects can safely be Ignored for thin blades. 


Comparison with Experimental Results 

In order to verify the present theoretical development of the equations 
of motion, a typical set of results were obtained for rotating, preconed, 
untwisted blade cases corresponding to those obtained from the NASA Lewis spin 
test rig (ref. 16). The test rig Is capable of accommodating bladed rotors up 
to 51 cm (20 In.) In diameter, which can be spun to 16 000 rpm. A rotor cap- 
able of holding two blades with adjustable sweep (precone) and setting angles, 
was designed for this spin rig. Tests were conducted for nonrotating, steel, 
flat plates having an aspect ratio, L/b, of 3 and thickness ratio, d/b, of 
0.05. Frequencies of the lowest five modes as given by the static spin rig 
tests were 94.2, 572, 586, (1643.5+11.5) and 1795 Hz respectively. Results 
obtained from the holography tests were respectively 93.0, 572, 586, 1628, and 
1793 for the lowest five mode frequencies. Of these frequencies, the first two 
correspond to the fundamental bending and torsion modes of the flat plate 
respectively. 

In order to match these frequencies from theoretical considerations, the 
elastic modulll were determined from standard relations. A value of 0.283 
lb/ln^ was chosen for the mass density. The effective root distance from 
the center of rotation of the spin rig varies slightly with the precone. For 
precones of zero, 22.5° and 45°, the effective root distances were 4.0, 3.952, 
and 3.816 In., respectively. Using these geometric and physical properties, 
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numerical results were generated for various setting angles, precones, and 
rotational speeds. Corresponding results obtained from the test rig are com- 
pared to the present theoretical results In table III. 

An examination of these results, presented In table III, Indicates that 
there Is a reasonable agreement between the theoretical and experimental 
results for all the cases considered. Furthermore, the trends shown by the 
theoretical and experimental results are consistent. From the comparison of 
the theoretical and experimental results shown In table III (of the order of 
+6 percent difference In most of the cases considered), and from the close 
agreement between the theoretical results and NASTRAN generated results pre- 
sented In tables I and II, It Is concluded that the present theoretical form- 
ulation Is accurate, and that It Is capable of producing accurate results for 
the parameters considered here. 


Vibration and Stability of Preconed, Rotating Blades 

In order to determine the Individual and combined Influence of pretwist, 
precone, Coriolis forces, second degree geometric nonl Inearltles and rotation 
on the coupled frequencies or stability boundaries, parametrlcs studies were 
conducted for various thickness ratio blade cases. A typical set of such 
results are shown In table IV. Also Included In this table are the results 
produced by MSC/NASTRAN. Figure 2 shows a further comparison of frequency 
ratios obtained from the present beam theory and from MSC/NASTRAN for typical 
precones and thickness ratio cases. A comparison of these results Indicates 
that there Is a good agreement between the present theoretical results and 
MSC/NASTRAN results for blade cases having moderate thickness ratios (d/b > 
0.10), moderate precones (Bgc = 15°) and for a wide range of rotational 
speeds (Q/o>-| up to 2.0). This trend of close agreement between the beam 
theory results Including second degree geometric nonllnearltles and those from 
MSC/NASTRAN continues for blades with Increasing thickness ratios for a wider 
range of precone angles and rotational speeds (see fig. 2(b)). In the absence 
of precone, the present beam theory results agree very closely with MSC/NASTRAN 
results for all thickness ratios, (see fig. 2(a)), over a wide range of rota- 
tional speeds. However, when the thickness ratio Is decreased, the agreement 
between the beam theory results and MSC NASTRAN results Is found to be close 
for low rotational speeds (up to about Q/u-j = 0.5 to 0.8) only (see fig. 2(c) 
and table IV). 

From the parametric results generated, It was found that the present beam 
theory predicts torsional Instability for thin blades (d/b = 0.05 or 0.06) at 
lower rotational speeds than does MSC/NASTRAN. This trend can be verified from 
the results presented In table IV for the case of a thin blade (d/b = 0.05) 
having 30° pretwist and 15° precone. For this particular blade configuration, 
the beam theory results are In close agreement with those of MSC/NASTRAN for a 
rotational speed parameter value of n/w-) = 0.5. However, when the rota- 
tional speed Is Increased, the torsional frequency predicted by the beam theory 
was found to be unstable, although the corresponding one from the MSC/NASTRAN 
calculation was stable. Similar trends were observed for this blade case even 
In the absence of pretwist (see fig. 2(c)). Furthermore, a comparison of the 
results presented In table IV for the present blade case with fi/u] = 1.0 
Indicates that the flexural frequencies predicted by the beam theory and 
NASTRAN calculations are still In very close agreement. Similar trends of 
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Instability could also be found for blades of the same thickness ratio, rota- 
tional speed and pretwist but larger values of precones. 

From the comparison of the results presented In table IV, the following 
observations are made: 

(1) For blades of moderate thickness ratios (d/b of the order of 0.10 
and greater), for moderate precones ( (3p C of the order of up to 15°), and 
for rotational speeds of up to fi/w| = 2.0, Inclusion of geometric non- 
linearities up through second degree only appears to be adequate for a fair 
prediction of the blade frequencies. Thus, the present second degree equa- 
tions appear to be adequate for application to helicopter blade vibration 
analysis. 

(2) The effect of Coriolis forces on the frequencies of rotating preconed 
blades Is Insignificant for low thickness ratio blade cases, but could become 
significant for high thickness ratio blade cases (refer to the results pre- 
sented In table IV). Conclusions concerning the linear and nonlinear Coriolis 
forces on the frequencies of rotating blades that were presented In refer- 
ence 11 are valid here also. 

(3) Thin blades possessing large precones and subjected to considerable 
rotational speeds exhibit torsional divergence at a much earlier rotational 
speed when the present beam theory Is used than Is the case when MSC/NASTRAN 
Is used. 

In order to acquire further Insight Into the torsional Instability of thin 
blades, results produced by the present beam theory and MSC/NASTRAN are pre- 
sented In table V for an untwisted, 15° preconed, thin blade (d/b = 0.05) for 
various rotational speeds. To facilitate a clearer understanding, the fre- 
quencies are listed In a columnwise fashion, each column corresponding to one 
of the bending modes (FI, F2, F3, F4) In the flatwise direction, the torsion 
mode (Tl) or the edgewise bending mode (SI) respectively. Also Included are 
the dimensionless steady state tip deflections (w,V, 4 >) In the last three 
columns of this table. An examination of the results presented In table V 
Indicates that for all rotational speeds considered, the flatwise steady state 
deflection, w, Is the most significant one while the edgewise and torsional 
deformations are Insignificant, as expected for the untwisted blade. Further- 
more, the steady state deflections produced by the present beam theory are In 
excellent agreement with those from MSC/NASTRAN. Since the flatwise deflection 
Is the only significant one that can contribute to the nonlinear terms In the 
perturbation solution, corresponding terms In the torsional equations causing 
Instability could easily be Identified. For the system to be unstable, the 
off-diagonal terms In the stiffness matrix should become large and negative. 

An examination of the stiffness matrix, [JK], for the perturbation solution, 
presented In appendix A Indicates that the^ nonlinear terms TS^j^, T9^jk 
and TlO^k are associated with the steady state equilibrium coordinates 
Wok, anc * ™ a t these terms arise from the terms discussed by Mil et al. 

(ref. 14). To verify this, results for same blade configuration presented In 
table V are reproduced In table VI, but with the Mil's terms discarded. From 
table VI, It can be seen that the Instabilities predicted by the present beam 
theory, shown In table V, vanish once the Mil's terms are discarded. However, 
the coupling trend between the fundamental torsional mode (Tl) and the funda- 
mental edgewise mode (SI), which now becomes evident by a mutual comparison of 
the results presented In tables V and VI [which decreases the lower uncoupled 
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frequency and Increases the higher uncoupled frequency with Increasing rota- 
tional speeds] Is also absent. Further, the present beam theory frequencies 
are more closely In agreement with those from MSC/NASTRAN when the Mil's terms 
are present than when they are absent, In so far as stable configurations are 
compared. Thus, It can be concluded that the Mil's terms are necessary for a 
reasonable prediction of the coupled bending-torsional frequencies, however, 
they also can cause torsional divergence for thin blades at larger precones and 
rotational speeds. It appears, then, that the level of approximation used for 
the Mil's terms may not be adequate, and Inclusion of the complete set of Mil's 
terms, (presently truncated In accordance with the ordering scheme or degree 
of geometric nonlinearity retained In the equations), may Improve the present 
theoretical prediction of Instability. These aspects are yet to be 
Investigated. 

Finally, the effect of various nonlinear terms on the coupled bendlng- 
bendlng-torslon frequencies of a particular blade with a 15° or 45° precone are 
presented In table VII. For assessing the Individual effect of each key non- 
linear terms, solution of the present nonlinear equations Is accomplished by 
retaining all nonlinear terms other than the one key element being addressed. 
Thus, In table VII, results under the column with T5^j|< = T6^|< = 0 represent 
nonlinear frequencies obtained by discarding the nonlinear terms arising from 
the tension effects In the torsion equation, results under the column with 
T1 1Jk = T4 1Jk = T0 1Jk = 0 correspond to elimination of those nonlinear terms 
associated with GJ In the two bending equations together with GJ(v'w")' In 
the torsion equation, and those under the column T3^j|< = T9^j|< = 0 (T2^j|< = 
TlO^jk = 0 for a = 0) correspond to elimination of the Mil's terms. Further, 
when A^jk Is set to zero, the foreshortening effects are Ignored, while 
D 1jk = E 1jk = 0 neglects the nonlinear terms arising from (Tw')' and (Tv 1 )'. 
From the results presented In table VII, and by comparing them with those 
presented In table IV one can observe that: 

(1) The affect of the terms (<t>'v'GJ)" and (<j>"w"GJ)‘ In the two bending 
equations (1) and (2) and the corresponding term (GJv'w")' In the torsion equa- 
tion (3) Is negligible on the coupled bendlng-bendlng-torslonal frequencies. 
These terms can therefore be neglected. 

(2) The affect of nonlinear terms arising from the tension coupling In 
the torsion equation Is also negligible on the coupled frequencies. 

(3) The Mil's terms are extremely Important In bringing In the bending- 
torsional coupling, and can also cause torsional divergence for thin blades. 
Further study with a careful consideration of these terms In their complete 
form Is necessary. 

(4) Compared to the Mil's terms just discussed, all other nonlinear terms 
that result due to the torsional coupling In the bending-torsion equations 
appear to be unimportant. 

(5) Effects of the nonlinear terms arising due to foreshortening, tension 
coupling with flexural motions, [(Tv')' and (Tw')'], and Coriolis effects were 
discussed In reference 11 In detail. These terms are extremely Important and 
must be retained In the equations. Conclusions concerning these terms as 
drawn In reference 11 are valid even In the presence of torsional coupling 
addressed In the present work. 
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CONCLUDING REMARKS 


The coupled bendlng-bendlng-torslon equations of dynamic motion of 
rotating linearly pretwisted and large preconed blades of symmetric cross 
section Including second degree geometric nonllnearltles and Coriolis effects 
are derived. These equations are solved by using the Galerkln method and a 
linear perturbation procedure. Natural frequencies and steady state deflec- 
tions produced by the solution of the present set of nonlinear equations are 
compared to those produced by MSC/NASTRAN calculations and also to those from 
experiment. Close agreement of the present theoretical results with those from 
other methods establishes the accuracy of the equations developed and the 
method of solution adopted. Parameter limits within which the second degree 
geometric nonllnearltles are adequate for a fair prediction of natural fre- 
quencies and steady state deflections are established. The following specific 
conclusions have emerged from the present Investigation. 

(1) For blades of moderate thickness ratios (d/b > 0.10) and moderate 
precones (B pc < 15°), Inclusion of geometric nonllnearltles up through sec- 
ond degree appears to be adequate for a fair prediction of the coupled fre- 
quencies and steady state deflections. 

(2) The present nonlinear equations Indicate torsional divergence at lower 
rotational speeds when the thickness ratio Is decreased and the blade precone 
Is Increased than do the finite element calculations. The nonlinear terms 
contributing to the coupling between bending and torsional motions, often 
referred to as kinematic pitch coupling terms and first discussed by Mil et al. 
(ref. 14) were found to be responsible for this torsional divergence. These 
terms are extremely Important In producing the accurate coupling between the 
bending and torsional frequencies. Inclusion of a complete set of these terms 
(presently appearing In the equations In a truncated form) may lead to a 
satisfactory stability calculations. 

(3) The affect of linear and nonlinear Coriolis forces on the coupled 
frequencies of thin blades Is found to be negligible. The Coriolis force terms 
can therefore be safely Ignored In analyzing advanced turboprop type blade 
configurations. However, the Coriolis effects must be retained In analyzing 
thick blades, as discussed In reference 11. 

(4) The affect of nonlinear terms arising from the tension coupling In the 
torsion equation, and the terms (GJ^'v 1 )", (GJ4>"w")' and (GJv'w")' appearing 
respectively In the flatwise, edgewise and torsion equations, Is found to be 
negligible on the coupled frequencies and steady state deflections for 
preconed rotating blades. 
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APPENDIX A: 


THE GALERKIN INTEGRALS AND MODAL EQUATIONS 

The various Integrals arising from the Galerkln process are defined 
below, and these are used In representing the modal equations In matrix forms 
subsequently: 

1 1 j 

5 1j = J dn = 5 e l e j dn = 1 j 1 = j 

C 1 /• T1 1 _ • _ _ 

A 1jk = J ♦i J dx dT i 


f l 

B ij = 3 dn 

; . 

c ij = J ^i^j 5 dn 


i „ i 


D 1jk = / ^l^j I x ) dx dn 

P « 

E 1 jk = J Vj'hc dn 

r 1 V / 2 

F ij = 3 *i*j y cos 0 


+ sin 2 e ) dn 
d^ 




1 „ 
= 3 vj 


sin 20 dn 


f " 

H. , = I q/ . t|r . cos 26 dn 
■J z ’ J 


1J 


i 

J 




J * ii i 

cos 20 dn 
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U 1 = J *1 sin 20 dr > 


T2 


Uk 


T3 


1 jk 


v. = f cos 29 dn 
1 o 1 

J / I I I I II II II* I \ 

T1 IJk = j *l(Vk + 2 Vk + e j *k) dn 

J* ’J' 1 j 0 j4' 1V sin 20 + 20j\jr k sin 20 + 0j<|> k sin 20 + 4Y0j'k k cos 20 

i ii 2 11 

+ 4Y0jt k cos 20 - 4 y 0j<l» k sin 20 J dn 

* J [v 


I II I 


cos 20 + 20jt)» k cos 20 + 0j »v k cos 20 - 4Y0j<l' k sin 20 


- 4Y0jl» k sin 20 - 4 y cos 20J dn 
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S8< = f e, sin 
1 o 1 


cos 2e dn 


0 cos 0 dn 
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The linear and nonlinear parts of the steady state equilibrium equations 
presented below in the matrix form, [JL + NL 3 f ^ 0 > =• <&> 
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The mass, gyroscopic and stifness matrices resulting from the pertur- 
bation equations are designated [M], [Cj and [K] respectively, and are 
presented in the following: 


- 

1 

« . . 1 0 
1J I 1 

r 

o 

0 1 s u 1 

_L_ 1J _J 

0 

o 

o 

] 

6 ij 

L 



o 

1 

2 sin ®Pc 5 ij - 

y 

~ 2 cos B Pc k w ok^ D ikj - E ikj) 

0 

- 2 sin 8 Pc 5 ij 

- 2 cos B Pc? w ok A ijk 

y 1 

- 2 c ° s Bp c ^'V ok (D ikj . - E ikj .) 

- 2 cos B Pc? v ok A ijk 

o 

L 

o 1 

0 

O 
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sin 2 8pc 4lJ - cos 2 • pc (B 1J - C. .) 


- sin 8 Pc cos 8 Pc k w ok A ijk + 


+ f ll ? q ’ok T2 ikj 


:?«ok[ D 


sin 8pc cos 8 Pc Y w ok Pijk + D ikj - E ijk - E ikj 


lKGi^rHij) 


sin 8 pc cos 6 pc ^v ok (0, kj - E ikj ) 


{? ! ij 


2 ^ J ij * 2y K iJ 


+ f 10? , ’ok T4 ikj 
+ f ll£*ok T3 1kJ 




sin 6 pc cos Bp c k v ok (T5 ikj " T6 ikj ) 

- sin 6pc cos e pc tOj, ♦ f 7 ? v o k T8 ikj 

" 2yE 8 s1n 8 Pc cos 8 Pc S6 ij 

- 2fj sin 6pc cos 8pc S2 1;j 

+ f 9? v ok T9 ikj - 2f 9? w ok T10 ikj 


- sin 8 Pc cos 8 Pc k v ok A ijk 


! i j + 

f io? 

f li?*ok T3 ikj 


2tCJ 1 j - YK ij 


f 2Z ti 

T 10 k *ok IA ikj 


j, ](^ ■ ‘) 


- s i j - cos (^(B^ - C 1 j) 


a Pc k w ok< D ijk “ E ijk) 


p Pc 


M ij - H? - V (G ij + YH ij } 


- f ^ ® T2 
T 11 k *ok Ikj 


f E w 

"7 k w . 


k~ w ok T8 ijk 


f l s1n a Pc cos a Pc S7 ij 


+ f 9? W 0k T9 ijk 
- 2f 9? V ok T10 ijk 


" f 10? v ok T1 1jk 
+ f ll^ w ok T2 ijk 


+ f li? v ok T3 ijk 


f 10? w ok T4 ijk 

V 

+ f ll k w ok T3 ijk 


_ f ll£ v ok T2 ijk 


f l cos 8 Pc S1 ij 
+ f 4 S4 ij + sin 8 p c cos 8 Pc 
£ w ok (T5 ijk " T6 ijk ) 


-cos e pc S5 ij 
+ 2Yf 5?W T7 ijk + T7 ikj) 

- R ij^ 2f 5 - f 6 * f 7 + f 8 cos2 
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APPENDIX B - NOMENCLATURE 


A i jk* B i y L v etc ' 

b. d 
d/b 


{BI, {X}, {X o I 


cross-sectional area of blade 
modal Integrals (see appendix A) 

breadth (chord) and thickness of blade 

thickness ratio 

vectors 


[£] 

E 

f r . • • f 

G 

1. J. k 

I 


11 


1 • 
nn 5? 


J 


kA 


km 

kml » km2 
L 


EL]. EJLN] 


m 

t"] 


n 

P 

R 

T 


modal damping matrix (gyroscopic matrix) 

Young's modulus 
coefficients (see eq. (13) 

shear modulus 
dummy Indices 

polar second moment of area about centroid 

area moments of Inertia about major and minor principal 
centroldal axes 

torsional stiffness constant 

blade cross section polar radius of gyration 

blade cross section mass radius of gyration 

principal mass radii of gyration 

length of beam 

Linear and nonlinear components of the matrix 
representing steady state equilibrium equations 

mass of blade per unit length 

modal mass matrix 

number of nonrotating modes for each of the flatwise 
bending, edgewise bending, and torsional deflections 

natural radian frequency 

radius of disc 

blade tension 
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t 

U, V, w 


time 

displacements of the elastic axis In X, Y, Z directions, 
respectively 


v, w, 9 

v ,, w ,, 9 . 
o j oj j 

X 

y, z 

a 

a r V 

Bpc 

y 

4v r 4 V **j 
4 u 

n 

0 

0j(n) 

e pt 

X I 
5 

p 

T 

^j(n) 

“1 

Q 

( )' 

C) 


dimensionless bending and torsional deflections 
steady-state equilibrium quantities 

running coordinate along X-axis 
centroldal principal axes of beam cross section 
setting angle (collective pitch) 
constants for assumed mode shapes 

precone angle 

total pretwist of the blade over Its length 
perturbation quantities 

Kronecker delta 

nondlmenslonal length coordinate, x/L 
geometric pitch angle, a + yn 
nonrotating torsional mode shape 
pretwist at a distance n from root, yn 
frequency parameter, V EI nn^AL 4 
nondlmenslonal rotational parameter, EI nn /pAL 4 Q 2 
mass density of blade material 
dimensionless time, nt 

nonrotating flatwise and edgewise bending mode shapes 

exact fundamental mode frequency of straight, nonrotating 
beam, 3.51602 V| 

rotor blade angular velocity, rad/sec 

primes denote differentiation with respect to x or n 

dot over a parameter represents differentiation with 
respect to t or t 
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TABLE I. - CONVERGENCE PATTERN OF FREQUENCY RATIOS (p /Xj_) OF A PRETWISTED, PRECONED, ROTATING BLADE 
INCLUDING SECOND DEGREE GEOMETRIC NONLINEARITIES AND CORIOLIS EFFECTS 


[n/o)l = 0.5, Bpc = 15’, Y = 30’, a = 0°, R = 0, d/b = 0.05, L/d = 200.] 


Mode 

number 

Results from perturbation solution: Galerkin method 

MSC NASTRAN 
(500 CQUAD4 
elements) 

n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

n = 6 

n = 7 

n = 8 

n = 9 

1 

6.7647 

3.9948 

3.9770 

3.9715 

3.9710 

3.9840 

3.9755 

3.9699 

3.9586 

3.9923 

2 

64.5292 

38.4098 

20.4184 

20.1127 

20.0309 

20.0334 

19.9880 

19.9823 

19.9707 

20.0707 

3 

73.7530 

55.3224 

52.5086 

52.3541 

52.4005 

53.3161 

51.8719 

51.6155 

51.4160 

57.7067 

4 


87.9097 

83.8938 

61.4177 

60.3757 

58.1215 

59.7126 

60.1018 

60.1548 

64.5911 

5 


206.4990 

137.3498 

85.8415 

85.3058 

84.9539 

86.2480 

86.5334 

86.7023 

86.1357 

6 

— 

438.1967 

205.3006 

205.2917 

123.9113 

123.4958 

121.5896 

120.8147 

120.6553 

121.1038 

7 

— 


345.6017 

274.0265 

205.4352 

198.2435 

198.2884 

199.2196 

199.0793 

199.4268 

8 


— 

457.0644 

343.8727 

343.0762 

343.5164 

205.8776 

204.6820 

204.3840 

219.5793 


TABLE II. - CONVERGENCE PATTERN OF STEADY-STATE TIP DEFLECTIONS 


[G/w-j = 0.5, 3p C = 15°, y = 30°, a = 0°, R = 0, d/b = 0.05.] 


Method 

Number of assumed 
modes or CQUAD4 
elements 

Steady-state tip deflection 
(center line deflection) 

w 

<1 

<P 

Galerkin 

n = 1 

-0.018471 

0.0018304 

-0.0044147 


n = 2 

-0.054886 

0.0065991 

-0.0016439 


n = 3 

-0.054606 

0.0058486 

-0.0012554 


n = 4 

-0.054841 

0.0058761 

-0.0011558 


n = 5 

-0.054832 

0.0058467 

-0.0011435 

MSC NASTRAN 

500 elements 

-0.054309 

0.0058372 

-0.0010209 











TABLE III. - COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


[d/b = 0.05, L/d = 60, y = 0°, L = 152.4 mm, E = 32.3xl0 6 psl, 
G = 13.8xl0 6 psl, p = 0.283 lb/1n 3 .] 


■ 

a, 

degree 

ft, 

rpm 

Method 


Frequency, 

Hz 


Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

45 

90 

1800 

(a) 

94.6 

573.4 

594.6 

1656.2 

1883.1 




(b) 

95.5 

571 .0 

599.0 

1657.0 

1798.0 



4200 

(a) 

95.7 

565.3 

620.3 

1669.7 

1824.4 




(b) 

100.0 

545.0 

640.0 

1650.0 

1823.0 



5400 

(a) 

95.8 

557.3 

641 .1 

1675.2 

1838.6 




(b) 

100.0 

522.0 

675.0 

1642.0 

1854.0 

0 

30 

3600 

(a) 


582.2 

623.9 

1683.4 

1825.8 




(b) 

129.0 

570.0 

630.0 

1688.0 

1807.0 



4800 

(a) 

146.5 

588.7 

647.0 

1700.4 

1837.4 




(b) 

146.0 

582.0 

665.0 

1716.0 

1817.0 



6000 

(a) 

168.6 

597.0 

673.4 

1714.1 

1853.3 




(b) 

166.0 

587.0 

706.0 



22.5 

90 

2400 

(a) 

Bn 

574.3 

603.5 

1665.4 

1815.8 




(b) 

IwEI 

571 .4 

602.9 

1664.0 

1799.0 



4800 

(a) 

■ 

573.3 

644.2 

1700.6 

1836.1 




(b) 

121 .2 

562.0 

654.0 

1684.0 

1826.0 



6000 

(a) 

131 .8 

572.3 

673.3 

1721 .9 

1855.1 




(b) 

134.5 

553.0 

688.0 

1695.0 

1854.0 

45 



(a) 

94.8 

570.2 

591 .0 

1653.1 

1780.5 




(b) 

95.9 

574.0 

588.0 

1655.0 

1799.0 




(a) 

96.3 

WEBm 

593.0 

1655.0 

1675.0 




(b) 

98.6 

Eli 

589.0 

1662.0 

1803.0 

22.5 



mm 


572.5 

593.4 

1655.3 

1796.3 




Imp 


574.0 



1798.0 




jW 

111.5 

489.9 

625.2 

1569.0 

1684.5 




1EM 

124.4 

561 .0 

626.0 




0 

0 

n 

nl|fl 

ESI 

mi 

1652.0 

1810.0 




wM 

■ 

Essl 

EH 

1655.0 

1795.0 


a Present theoretical results. 
^Experimental result. 















































TABLE IV. - COMPARISON OF FREQUENCY RATIOS OF PRETWISTED, PRECONED, ROTATING BLADES OF 

VARIOUS THICKNESS RATIOS 


[a = 0°, R = 0, L/b = 10, geometric nonllnearltles Included.] 


jjj^^ 

Y t 

degree 

0p C . 

degree 

Q/oy| 



Method Frequency ratio, 

p/V | 



Mode 2 


Mode 4 

Mode 5 




15 

0.8 

(a) 


23.1186 

30.4598 

62.7814 







(c) 


23.1970 

31 .9459 

62.8686 


122.3435 




1.0 

(a) 

5.1488 

23.7212 


63.3946 

68.1157 

122.6629 





(c) 


23.7841 

30.7285 

63.3849 

73.5597 

122.8331 




2.0 

(a) 

8.3781 

21 .3722 

28.3088 

68.3325 


127.9158 





(c) 

8.3327 

28.3285 

28.7811 


75.6988 


0.25 

0 

45 

0.8 

(a) 

3.9975 


22.7272 

61 .8443 

62.4559 






(b) 






mSmmm 





(c) 

4.2598 



61 .3970 

67.2376 

89.0158 

0.25 

0 

45 

1.0 

(a) 

4.4286 


23.2294 


62.9189 






(b) 

4.6184 



59.4123 

62.8856 

92.1343 





(c) 

4.8122 




67.1483 

91 .2117 

0.06 

0 

15 

0.5 

(a) 

3.9654 

22.4572 

48.9361 

62.1175 

72.6927 

121.3362 





(c) 

3.9866 

22.5073 

52.2357 

62.1740 

75.2841 

121 .4498 

0.06 

0 

45 

0.5 

(a) 

3.6499 

22.2620 

Unstable 

61.9376 

75.9085 

121.1543 





(c) 


22.2822 

40.7961 

61 .7824 

85.4823 

121 .0251 

0.05 

30 

15 

0.5 

(a) 

3.9710 

20.0309 


60.3757 

85.3085 

123.9113 





(c) 

3.9923 

20.0707 

57.7067 

64.5911 

86.1357 

121 .1038 



15 

1.0 

(a) 

5.1659 

21.7959 

Unstable 


86.4133 

125.2654 





(c) 

5.1656 



61.8346 

92.8652 

121 .8441 


30 

15 

1.0 

(b) 

5.1411 

15.4978 



67.1784 

113.9297 





(c) 

5.1619 

15.5424 

26.0745 

60.7458 

70.4579 

109.2493 

0.20 


45 

0.8 

(b) 


12.9163 

23.3475 

58.2695 

64.4610 

114.5998 





(c) 

4.2651 

14.2839 

24.4469 

58.7436 

73.0801 

110.8504 

0.20 

30 

45 

1.0 

(b) 

0.1533 

7.2997 


54.6244 

63.4657 






(c) 

Unstable 

9.1287 

18.2334 


84.9760 


0.05 

30 

45 

0.8 

(a) 

Unstable 










(c) 

Unstable 







a Present beam theory Including Coriolis effects 
b Present beam theory Ignoring Coriolis effects. 
c MSC NASTRAN results Ignoring Coriolis effects. 













TABLE V. - COMPARISON OF NONLINEAR FREQUENCY RATIOS AND STEADY-STATE DEFLECTIONS OF A PRECONED BLADE AT VARIOUS 

ROTATIONAL SPEEDS 


[a = 0*, y = 0°, R = 0*, 8 pc = 15* , L/b = 10, d/b = 0.05, Coriolis effects not included. Unstable as negative 
p 2 obtained for I torsional frequency.] 


n/ u>i 

Method 

Frequency ratio, p/xi 

Steady 



Mode 1 (FI) 

Mode 2 (F2 ) 

Mode 3 (F3) 

Mode 4 (Tl) 

Mode 5 (SI) 

Mode 6 (F4) 

w 

0.1 

(a) 

3.5344 

22.0513 

61.7139 

68.5155 

70.5703 

120.9191 

-0.002803 


(b) 

3.5594 

22.1618 

61.9851 

68.1727 

70.1475 

121.4925 

-0.002767 

0.3 

(a) 

3.6805 

22.1860 

61.8476 

62.9968 

74.5192 

121.0572 

-0.02326 


(b) 

3.7062 

22.2934 

62.1044 

63.9497 

74.0049 

121.6067 

-0.02297 

0.5 

(a) 

3.9654 

22.4573 

62.1175 

50.3754 

79.1494 

121.3362 

-0.05572 


(b) 

3.9950 

22.5585 

62.3354 

57.2115 

79.8750 

121.8374 

-0.05500 

0.8 

(a) 

4.6098 

23.1186 

62.7814 

Unstable 

82.3113 

122.0243 

-0.10620 


(b) 

4.6451 

23.2026 

62.8837 

48.1742 

86.2660 

122.3773 

-0.10463 

1.0 

(a) 

5.1448 

23.7213 

63.3947 

Unstable 

82.2947 

122.6629 

-0.13383 


(b) 

5.1781 

23.7926 

63.4078 

44.2514 

88.2218 

122.8856 

-0.13173 

2.0 

(a) 

8.3781 

28.3089 

68.3326 

Unstable 

81.0269 

127.9159 

-0.20305 


(b) 

8.3385 

28.3445 

68.0883 

39.0512 

87.6030 

127.5273 

-0.20001 

3.0 

(a) 

11.8998 

34.6488 

75.8305 

Unstable 

82.1026 

136.2440 

-0.22464 


(b) 

11.7580 

34.6610 

75.5539 

34.6610 

84.4642 

135.4860 

-0.22205 


Steady-state tip deflections 


a Present beam theory. 

b M$C NASTRAN results ignoring Coriolis effects. 


TABLE VI. - FREQUENCY RATIOS AND STEADY-STATE DEFLECTIONS OF A PRECONED BLADE 

NEGLECTING MIL'S TERMS 


[a = y = 0°, B pc = 15°, R = 0, d/b = 0.05, L/b = 10, T2ij k = T3ij k = 
T9lj k = T10 1jk = 0- Coriolis effects ignored.] 


Q/t o*| 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

Mode 6 

w 

H 

<P 

0.1 

3.5344 

22.0513 

61 .7139 

68.7889 

70.3204 

120.9191 

-0.002803 

_ 

_ 

0.3 

3.6805 

22.1860 

61.8476 

68.8022 

70.3218 

121 .0572 

-0.02326 

- 

- 

0.5 

3.9654 

22.4573 

62.1175 

68.8274 

70.3263 

121 .3362 

-0.05572 

- 

- 

0.8 

4.6098 

23.1186 

62.7814 

68.8879 

70.3394 

122.0243 

-0.10620 

- 

- 

1 .0 

5.1448 

23.7213 

63.3947 

68.9454 

70.3505 

122.6629 

-0.13383 

- 

- 

2.0 

8.3781 

28.3089 

68.3326 

69.4571 

70.4161 

127.9159 

-0.20305 

- 

- 

3.0 

11.8998 

34.6488 

70.2570 


75.8305 

136.2440 

-0.22464 

- 

- 




















TABLE VII. - EFFECT OF IGNORING CERTAIN NONLINEAR TERMS ON THE FREQUENCY RATIOS, p/x lf OF ROTATING, PRECONED BLADES 


[d/b = 0.06, L/b = 10, a = y = 0* , R = 0, n/w \ = 0.5. Coriolis effects included.] 


V* 

degree 

Mode 

number 

T5 i jk " T6 ijk “ 0 

T1 ijk " T4 1 jk ' 
T8 1jk ■ 0 

T3 i jk " T9 ijk " 0 

Tl^ through 
T10 ijk are 0 

A uk - °* T1 uk 
through TIO^ 

are 0 

A 1jk jk = ^ i jk 8 

T1 1jk t0 T10 ijk ■ 0 

All nonlinear 
terms ignored 

15 

1 

3.9650 

3.9650 

3.9651 

3.9651 

3.9592 

3.9475 

3.9475 


2 

22.4572 

22.4572 

22.4572 

22.4572 

22.4563 

22.4497 

22.4497 


3 

48.9396 

48.8668 

58.6064 

58.6070 

58.6092 

58.6104 

58.6104 


4 

62.1176 

62.1176 

62.1186 

62.1186 

62.1171 

62.1113 

62.1113 


5 

72.6937 

72.7639 

68.8350 

68.8341 

68.8342 

68.8342 

68.8342 


6 

121.3362 

121.3362 

121.3362 

121.3362 

121.3360 

121.2972 

121.2972 

45 

1 

3.6558 

3.6557 

3.6474 

3.6474 

3.6171 

3.5543 

3.5543 

! 

2 

22.2614 

22.2614 

22.2614 

22.2614 

22.2573 i 

22.2251 

22.2251 


3 

Unstable 

Unstable 

58.6235 

58.6273 

58.6368 

58.6426 

58.6426 


4 

61.9371 

61.9371 

61.9412 

61.9412 

61.9360 

61.9076 

61.9076 


5 

75.9208 

76.1321 

68.8184 

68.8141 

68.8141 

68.8141 

68.8141 


6 

121.1544 

121.1544 

121.1544 

121.1190 

121.1187 

121.1255 

121.1255 








Beam theory 
divergence 



.2 .4 .6 .8 

Rotational speed parameter, Q/u 1 


[c)-jj ■ 0.5, a * y* 0°, (3 pc " 15°, R- 0 
Comparison of frequency ratio 
Figure 2. - Concluded. 
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